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Abstract 

' We present an algorithm for isolating the roots of an arbitrary complex polynomial p that also works 

for polynomials with multiple roots provided that the number k of distinct roots is given as part of 
the input. It outputs k pairwise disjoint disks each containing one of the distinct roots of p, and its 
multiplicity. The algorithm uses approximate factorization as a subroutine. 

In addition, we apply the new root isolation algorithm to a recent algorithm for computing the topol- 
(/3 ' ogy of a real planar algebraic curve specified as the zero set of a bivariate integer polynomial and for 

O I isolating the real solutions of a bivariate polynomial system. For input polynomials of degree n and bit- 

size T, we improve the currently best running time from d{r?z + r?x-) (deterministic) to 0(rfi + n^'c) 
(randomized) for topology computation and from 0(n^ +n^T) (deterministic) to 0{n^ +/i^t) (random- 
^ ized) for solving bivariate systems. 

o 

00 ! 1 Introduction 

In the main part of this paper, we give an algorithm for isolating the roots of a polynomial p{x) = L"=o P<^' 

■ with arbitrary complex coefficients that works for polynomials with multiple roots, provided that the num- 
I ber k of distinct roots is given as part of the input. More precisely, let zi, Zk be the distinct roots of p, 

let rrij := mult{zi,p) be the multiplicity of z;, and let :^ a{zi,p) '■= miny^,- — Z; | be the separation of Zi 
from the other roots of p. Then, the algorithm outputs isolating disks Ai = A{z\,R\) to = A{zk,Rk) for 
the roots of p as well as the corresponding multiplicities mi to m^. The radii satisfy /?, < ^ for all /, thus 
I the center zi of A, approximates z, to an error of less than If the number of distinct roots of p differs 

■ from k, we make no claims about termination and output. 

The coefficients of p are provided by oracles. On input L, such an oracle essentially returns binaiy 
fraction approximations /5/ of the coefficients pi such that \\p — Y!i=QPix'\\ ^ \\p\\. Here, \\p\\ := \\p\\ ^ = 
\pq\ + ■■■ + \p„\ denotes the one-norm of p. The details are given in Section lZTI 

The algorithm has a simple structure. We first use any algorithm (e.g. |l4l[l8][l6l|2T|) for approximately 
factorizing the input polynomial. It is required that it can be run with different levels of precision and that, 
for any given integer b, it returns approximations zi to z„ such that 



X 



<2-'\\p\\- (1) 



In a second step, we partition the root approximations z\ to z„ into k clusters Ci , . . . ,Q based on geometric 
vicinity. We enclose each cluster C, in a disk D, — A(z, , r, ) and make sure that the disks are pairwise disjoint 
and that the radii r, are not "too small" compared to the pairwise distances of the centers f,Q In a third 
step, we verify that the «-times enlarged disks A, = A(|/,/?,) = A(z/,« • r,) are disjoint and that each of 
them contains exactly the same number of approximations as roots of p counted with multiplicity]! If the 
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'This will turn out to be crucial to control the cost for the final verification step. For details, we refer to Sections l2.2.2l and l2.2.3l 
^It seems to be artificial to consider enlarged disks A, instead of the disks D, . This is mainly due to the usage of Rouche's Theorem 
to verify roots. 
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clustering and the verification succeed, we return the disks Aj , . . . , A^. and the number of approximations z G 
{zi , . . . ,£„} in the disk as the multiplicity of the root isolated by the disk. If either clustering or verification 
does not succeed, we repeat with a higher precision. 

If Pan's algorithm llT6l is used for the approximate factorization step, then the overall algorithm has bit 
complexitjH 

O ^n' + «2 £ logMfe) + n log (M((jr'"' ) • MiPr ' )) j (2) 

where Pj := Ylj^iizi — Zj)"'j = ^ and M{x) :~ max(l, |x|). Observe that our algorithm is adaptive 

in a very strong sense, namely, the above bound directly depends on the actual multiplicities and the 
geometry (i.e. the actual modulus of the roots and their distances to each other) of the roots. There is also 
no dependency on the size or the type (i.e. whether they are rational, algebraic or transcendental) of the 
coefficients of p. 

Our algorithm can also be used to further refine the isolating disks to a size of 2^'^ or less, where K is 
a given integer The bit complexity for the refinement is given by the bound in (|2]) plus an additional term 
0{n-K- max,- m,- ) . In particular for square-free polynomials, the amortized cost per root and bit of precision 
is one. 

For the benchmark problem of isolating aU roots of a polynomial p with integer coefficients of absolute 
value bounded by 2"^, the bound in ^ becomes d{ir' +n^z). The bound for the refinement becomes 
d{n^ + n^Z + nK), even if there exist multiple roots. 

For a square-free integer polynomial p, we are aware of only one method ||9] Theorem 3. 1] that achieves 
a comparable complexity bound for the benchmark problem. That is, based on the gap theorem from 
Mahler, one can compute a theoretical worst case bound bo of size 0(«t) with the property that if « points 
Zj e C fulfill the inequality ([T]) for a b> bo, then they approximate the corresponding roots zj to an error 
less than (7//(2n); cf. Lemma|4]for an adaptive version. Hence, for b > bo, Pan's factorization algorithm 
also yields isolating disks for the roots of p using (5(«^t) bit operations. Note that this approach achieves a 
good worst case complexity, however, for the price of running the factorization algorithm with b = 0(«t), 
even if the roots are well conditioned. In contrast, our algorithm turns Pan's factorization algorithm into 
a highly adaptive method for isolating and approximating the roots of a general polynomial. Also, for 
general polynomials, there exist theoretical worst case bounds |18, Section 19] for the distance between 
the roots of p and corresponding approximations fulfilling ([T). They are optimal for roots of multiplicity 
n(n) but they constitute strong overestimations if all roots have considerably smaller multiplicities. For 
the task of root approximation, the complexity of our method adapts to the highest occurring multiplicity 
(with bit complexity (5(nmax,OT,- • k:) for K dominating), whereas this is not given for the currently best 
method jTSI (with complexity (9(«^K')). 

We would also like to remark that we are aware of only one previous root isolation algorithm ifTSi that 
can cope with multiple roots, however, with at most one multiple root. In addition, the number of distinct 
complex roots as well as the number of distinct real roots must be given as an additional input. 

Finally, we aim to stress the importance of our root isolation method by applying it to the problems 
of computing the topology (in terms of a cylindrical algebraic decomposition) of a real planar algebraic 
curve specified as the zero set of an integer polynomial and of isolating the real solutions of a bivariate 
polynomial system. Both problems are well-studied |[T][n][l2l[T9]|7]|5][3l[8l [H. More specifically, 
we apply our method to a recent algorithm |:3] for computing the topology of a planar algebraic curve. 
This yields bounds on the expected number of bit operations which improve the currently best (which are 
both deterministic) bounds ||8][T4l from 0(n^T + n^T^) to(5(n^+n^T) for topology computation and from 
d{n^ + I'P x) to (5(n^ + iv't) for solving bivariate systems. 

indicates that we omit logarithmic factors. 
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2 Root Approximation 
2.1 Setting and Basic Properties 

We consider a polynomial 

pix)=p„x" + ...+poeC[x] (3) 
of degree n>2, where p,, ^ 0. We fix the following notations: 

• M(;c) := max(l, |x|), forx e M, 

• Tp denotes the minimal non-negative integer with < for all ; = 0, ...,«— 1, 

• IIpII •= IIpII 1 ■= IpoI + • • • + \pn \ denotes the 1-norm of p, 

• zi , . . . G C are the distinct complex roots of p, with k <n, 

• nii := mult{zi,p) is the multiplicity of Zi, 

• Gi:=G{zi,p) := miny^,- |zi — Zjl is the separation of Zi, 

• := M(max, log |z, |) the logarithmic root bound of p, 

• Mea(/?) — \pn\ ■Y\iM{zi)"^'' the Mahler Measure of p. 
The quantities Tp, Yp, \p„\ and Mea(/?) are closely related. 

Lemma 1. Fp < 1 + Xp and Xp — n — l < log ^^^j^jj^ < "Fp. 

Proof. By Cauchy's root bound max; | < 1 +max; \pi \ / \p„\, and thus max,Tog |z, | < 1 + Tp. Since Tp > 0, 
by definition, we have Fp < 1 + Tp. The ;-th coefficient of p is smaller than or equal to (") Mea(;?) < 

2"Mea(;') < 2"^^p+^\ Thus, from the definition of Tp, either Tp = or 2"^^^ > max,- > 2^p-^ □ 

We assume the existence of an oracle which provides arbitrary good approximations of the polynomial 
p. Let L > 1 be an integer We call a polynomial p = /5„jc" + . . . + po, with pi — s,- • 2^^ and i,, £ e Z, an 
approximation of precision L of p if \pi — pi\ < 2^^^'°^^"+'' £ < L+ [log(« + 1)] — [log||p||J, and 
log <L+ [log(« + 1)] + 1 for all /. When considering pi as infinite bitstring pf = sgn(p;) • Lit!°_oo ^<r2'"', 
bii £ {0, 1}, then we can obtain pi from the partial string which starts at index = [log||p||J and ends 
at index k2 = [log \\p\\\ L [log(n + 1)] , that is, i,- := 2' • sgn(/:',) • Lj^La-t ^A:2*^, and I = L + [log(« + 
1)] — [log||/?|jJ. We assume that we can ask for an approximation of precision L of p at cost 0{n{L + 
log«)) = 0{tiL). This is the cost of reading the approximation of precision L. The next Lemma states some 
elementary properties of approximations of precision L. 

Lemma 2. Let p be an approximation of precision L of p and L > 1. 

• ||/5||/2<||p|| <2||/5||. 

• IfL > Tp +4, then 2-^-l°g("+l) < |/5„| /4. 

• //2-^-i°s("+') IIpII < |/5„|/4, then \p„\/2 < \p„\ <2\p„\. 

Proof. From the inequality |||/5|| — ||p|| | < — we conclude that — H/^H | <2^^j|;9|| < j|;9|l/2. 
This establishes the first claim. For the second claim, we observe that 

|/5„-p„| <2-^-'°g("+')||p|| <2-^+V|p„| < \p„\/l6. 

Thus, \p„\ > \p„\/2 and 

2-L-iog(«+i) ij^jj < 2.2-^-'°g("+') ||;;|| < \p„\/8 < |p„|/4. 
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The third claim follows from 

\pn-Pn\ < 2-^-'°g(«+') Ibll < 2.2-^-i°s("+') 11/511 < \p„\/2. 

□ 

Lemma|2]suggests an efficient method for estimating />„. We ask for approximations p of precision L 
of p for L= 1,2, 4,... until the inequality 2-^-'°g(''+i) ||p|| < |p„|/4 holds. Then, \p„\/2< \p„\ <2|/5„|by 
part 3 of the Lemma. Also L < 2(Tp +4) by part 2 of the above Lemma. The cost is d{nZp) = (5(«^rp) bit 
operations, where we used the upper bound for Tp from Lemma [T] Observe that this bound depends only 
on the geometry of the roots (i.e. the actual root bound Fp) and the degree but not (directly) on the size 
of the coefficients of p. We next show that a "good" integer approximation F of Fp can also be computed 
with (5(n"Fp) bit operations. 

Theorem 1. An integer F e N with 

Fp<F<81og« + Fp (4) 

can be computed with d{n^Tp) bit operations. The computation uses an approximation of precision L of p 
with L = (9(«Fp). 

Proof. We first compute an approximation p of precision L of p with |/5„| /2 < |p„| < 2 • \p„\ as described 
above. Let K := 2L'°gl/^«IJ-i. Then K < \p„\/2 < \p„\ and K > 2Liog|p«/2|J-i > \p„\ /S. Consider the scaled 
polynomial q := p/K. Its leading coefficient lies between 1 and 8. Moreover, F^ = Fp, = Tp, and an 
arbitrary approximation of precision L of p also yields an approximation of precision Lof q since division 
by X is just a shift by log K bits and ||^|| = ||/7|| / k". Hence, we may assume that 1 < Ifnl < 8. 
Consider the Cauchy polynomial 

p{x):^\pnW~"Y^\pi\x' 

i=0 

of p. Then, according to Proposition 2.51], p has a unique positive real root E, e M+, and the following 
inequaUty holds: 

n 

max | < ^ < — — • max|z, | < 2n ■ max|z, | . 
I In 2 (■ ! 

It follows that p{x) > for all x > £, and p{x) < for all x < ^. Furthermore, since p coincides with 
its own Cauchy polynomial, each complex root of p has absolute value less than or equal to |^|. Let 
ko be the smallest non-negative integer k with /5(2*) > (which is equal to the smallest k with 2*' > ^). 
Our goal is to compute an integer F with ko <r < ko + I. Namely, if F fulfills the latter inequality, 
then M(maxi |z,|) < M{^) <2^ < AM{B,) < 8n-M(max; |z,|), and thus F fulfills inequality In order 
to compute a F with A:o < T < + 1, we use exponential and binary search (try k — 1,2,4,8, . . . until 
/7(2^) > and, then, perform binary search on the interval k/2 to k) and approximate evaluation of p at 
the points 2^: More precisely, we evaluate /5(2*^) using interval arithmetic with a precision p (using fixed 
point arithmetic) which guarantees that the width w of *B(/?(2*^),p) is smaller than 1, where f8(£',p) is the 
interval obtained by evaluating a polynomial expression E via interval arithmetic with precision p for the 
basic arithmetic operations; see ifTS] Section 4] for details. We use ifTSl Lemma 3] to estimate the cost for 
each such evaluation: Since p has coefficients of size less than 2'^p\p„ \ < 2^?+^, we have to choose p such 
that 

2-''+2(« + l)^2^''+^+"*-< 1 

in order to ensure that w < 1. Hence, p is bounded by 0{Tp + nk) and, thus, each interval evaluation 
needs (5(«(Tp + nk)) bit operations. We now use exponential plus binary search to find the smallest k such 
that *B(/5(2'^),p) contains only positive values. The following argument then shows that A:o < A: < A:o + 1: 
Obviously, we must have k> ko since p{2'') < and /5(2^) e *B(/7(2^),p) for all k < ko. Furthermore, the 
point jc = 2*0+^ has distance more than 1 to each of the roots of p, and thus |;5(2'''''+')| > |p„| > 1. Hence, 
it follows that 03(^(2*^" + l),p) contains only positive values. For the search, we need 

O(log^o) = O(loglog^) = 0(log(log« + Fp)) 
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iterations, and the cost for each of these iterations is bounded by d{n{Tp + nko)) = Oiir^Tp) bit operations. 

□ 

2.2 Algorithm 

We present an algorithm for isolating the roots of a polynomial p{x) = YJUoPi^' = /'nllf i (x — Zi)'"', where 
the coefficients pj are given as described in the previous section. We may assume A: > 1 ; the problem is 
trivial otherwise. If A: = 1, ~-pn-\/ {npn) is the root of multiplicity n. The algorithm uses some polynomial 
factorization algorithm to produce approximations for the roots zi , . . . ,Zjt, and then performs a clustering 
and certification step to verify that the candidates are of high enough quality. For concreteness, we pick 
Pan's factorization algorithm lITSl for the factorization step, which also currently offers the best worst case 
bit complexity. If the candidates do not pass the verification step, we reapply the factorization algorithm 
with a higher precision. Given a polynomial p with |z, | < 1 for \ <i<k, and a positive integer b denoting 
the desired precision, the factorization algorithm computes n root approximations £i , . . . ,£„. The quality of 
approximation and the bit complexity are as follows: 

Theorem 2 (Pan 11611 ). Suppose that \zi\ < 1 for I < i < k. For any positive integer b > «log«, complex 
numbers £i , . . . ,£„ can be computed such that they satisfy 

iiP-A,n"=i(-^-2')ii<2-iipii 

using 0{n) operations performed with the precision of 0(b) bits (or 0(bn) bit-operations). We write 
p := p„ ri/Li (x ~ Zi)- The algorithm returns the real and imaginary part of the Zi 's as dyadic fractions of 
the form A ■ 2^^ with A G Z, B e N and B = 0(b). All fractions have the same denominator 

The parameter b controls the quality of the resulting approximations. Note that Pan's algorithm requires 
all roots of the input polynomial to lie within the unit disk A(0, 1 ). Hence, in order to apply the above result 
to our input polynomial, we first scale p such that the roots come to lie in the unit disk. That is, we compute 
a r as in Theorem[T] and then consider the polynomial f(x) := p(s -x) = LlLo/i-*^' with s := 2^. Then, 
f(x) has roots = Zi/s G A(0, 1), and thus we can use Pan's Algorithm with b' nF + b to compute an 
approximate factorization /(x) := Y!i=ofix' •= fnl]'i=i(x — ii) such that ||/ — /|| < 2"''' ||/||. Let z,- := i- 1,- 
for all i and p(x) := p„ ■ UU - ^z) = /(^/^) = iLo//^'-^'- Then 

\\P p\\ = t \fi/'' /' A'l <t\fi'f-\< 2-'' t \M 

1=0 (=0 j=0 

<2-^"f |/,-A'|=2-^|1p||. 

1=0 

For the factorization of /, we need an approximation of precision b' of /, and thus an approximation of 
precision L of p with L = 0(nS + b) = d(nTp + b). The total cost is d(n^r + nb) bit operations. We 
summarize in: 

Corollary 1. For an arbitrary polynomial p and an integer b > «log«, complex numbers zi, - ■ ■ ,Zn can be 
computed such that 

\\p-p.^nu^-^^')\\<2-'\\p\\ 

using d(n^r + bn) bit-operations. We write p :~ PiiYl"=i(x ~ Zi). The algorithm returns the real and 
imaginary part of the Zi's as dyadic fractions of the form A ■ 2^^ with A G Z, B G N and B = 0(b + nFp). 
All fractions have the same denominator. 

We now examine how far the approximations zi, . . . ,£„ can deviate from the actual roots for a given 
value of b. Let A(z,r) be the disk with center z and radius r and let bdA(z,r) be its boundary. We further 
define P,- :== Y\j^i(zi - Z;)'"'. Then, p^"'-\zi) = mi\p„Pi- 

Lemma 3. Ifr < Oi/n, then 

\p(x)\>r''^^\pnPi\/4 

for all X on the boundary ofA(zi,r). 
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Proof. We have 



\p{x)\^\p„\\x~zirY\\x-zr 

> \pn\ i^-z,-i'"'ni^'-~^^'i"'' ■ (1 - v-mzi-z^r- 

> r'"' (1-1 /«)«-"" \p„ \l\\zi- Zj \"y > r"" \p„P, I /4. 



□ 



Based on the above Lemma, we can now use Rouche's theorem0EI to show that, for sufficiently large 
b, the disk A(z;,2^^/f^'"'^) contains exactly m,- root approximations. 

Lemma 4. Let p be such that ||p — /3|| < 2"* If 

b > max(8«,nlog(n)), andb is a power of two (5) 
2-^/(2"'') < l/(2«2), (6) 
2-b/{2m,) < (j./(2„), and (7) 

2-b/2 < \^ (gx 

- l6in + l)2^i'M{ziY 

for all i, the disk A(zi,2^^^^^'"'^) contains exactly m; root approximations. For i ^ j, let Zi and Zj be 
arbitrary approximations in the disks A(z,-,2^''/(^™'') and A{zj,2^^^^'^'"j')), respectively. Then, 

(l--)|zi-Z;| < \zi~Zj\ < {l + -)\zi-Zj\. 

n ' •'II •'I fi I •'I 

Proof. Let 

'l6 • (« + 1) • 2-''2'p \Pi\-^M{zi)") . 



It is easy to verify that 5; < 2^''/'^'"'' < min(l, (7,-)/(2n). The first inequality follows from (|8) and the 
second inequality follows from (|6]l and Q. We will show that A(z,-, 5,) contains m, approximations. To this 
end, is suffices to show that \{p — p)(x)\ < \p{x)\ for all x on the boundary of A(z,-, 5,). Then, Rouche's 
theorem guarantees that A(z,-, 5,) contains the same number of roots of p and p counted with multiplicity. 
Since Zj is of multiplicity m; and 5; < C7,-/«, the disk contains exactly m, roots of p counted with multiplicity. 
We have (note that |x| < (1 + 1 /(2n^)) ■ M{zi) for x e bdA(z,, 5,-)) 

\{p-p){x)\<\\p-p\\-\^T <2-'\\p\\\A" 

<2-'\\p\\-{\ + \/{2n^)r-M{zir 
<A-2-'-2'''\pn\-{n + \)-M{zir 
<5r\pnPi\/A<\p{x)l 

where the inequality in line three follows from ||p|| < (« + 1)|/9„|2^'', the first one in line four follows from 
the definition of 5,, and the last inequality follows from Lemma|3] It follows that A(z,,2^*''(^'"'') contains 
exactly m, approximations. Furthermore, since 5, < C7,/(2«) for all the disks A(z,, 5,), I <i <k, are 
pairwise disjoint. 

For the second claim, we observe that \zt —zt\ < 2^''/^^'"'' < a>/(2«) < \zi —Zj\/ (2«) for i = i,j and 
hence |£, — z,| + [fy — Zj| < Z;|/n. The claim now follows from the triangle inequality. □ 



Rouche's theorem states that if / and g are holomorphic functions with \ {f — g){^)\ < 1/(^)1 for all points .v on the boundary of 
some disk A, then / and g have the same number of zeros (counted with multiplicity) in A. 

'The reader may wonder why we are not using Gershgorin circles to cluster and veiify root approximations. The problem is that 
very good approximations of a multiple root generate very large Gershgorin circles which may then fail to isolate the approximations. 
Indeed, if two approximations are equal, the corresponding disks have infinite radius. 
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We have now established that the disks A(z,-,2-''/(2'«')), 1 < i < k, are pairwise disjoint and that the 
;-th disk contains exactly m, root approximations provided that b satisfies ^ to Unfortunately, the 
conditions on b are stated in terms of the quantities m,, (7, and | which we do know. Also, we do not 
know the center Zi- In the remainder of the section, we will show how to cluster root approximations and to 
certify them. We will need the following more stringent properties for the clustering and certification step. 

2-^/** < min(l/16,|P,-| /{{n + I) -2^"^"+^")) 

Let bo be the smallest integer satisfying (|5]l to ( fTOl i for all /. Then, 

b() = 0{nlogn + nrp + maXi{milogM{ar^)+\ogM{Pr^)). 

We next provide a high-level description of our algorithm to isolate the roots of p. The details of the 
clustering step and the certification step are then given in Sections [2.2.2| and |2.2.31 respectively. 

2.2.1 Overview of the Algorithm 

On input p and the number k of distinct roots, the algorithm outputs isolating disks A,- = A(z,-,7?,) for the 
roots of p as well as the corresponding multiplicities m,-. The radii satisfy Ri < Oi/ (64«). 

The algorithm uses the factorization step with an increasing precision until the result can be certified. 
If either the clustering step or the certification step fails, we simply double the precision. There are a 
couple of technical safeguards to ensure that we do not waste time on iterations with an insufficiently large 
precision (Steps 2, 5, and 6); also recall that we need to scale our initial polynomial. 

1. Compute the bound 2^ for the modulus of all roots of p, where T fulfills Inequality (|4|i. According 
to Theorem[T] this can be done with d{n^Tp) bit operations. 

2. Compute a 2-approximation X^2'^,lx G Z, of 1 1/9 1 1 / \p„ \ . According to Lemma|2]and the subsequent 
remarks, we can compute A with 0{n^rp) bit operations. 

3. Scale p, that is, f{x) := p{s-x), with s := 2^ , to ensure that the roots ^/ = Zi/S, i — 1, . . . ,k, of f are 
contained in the unit disk. Let b be the smallest integer satisfying (|5]l 

4. Run Pan's algorithm on input / with parameter b' :=b + nY to produce approximations , . . . , 4, for 
the roots of /. Then, £, s ■ t,i are approximations of the roots of p, and ||/3 — p\\ < 2^^* where 

P{x) ■■= Pnl\"=lix-Zi). 

5. If there is a £, with z, > 2""+^ return to Step|4]with b := 2b. 

6. IfU"=iM{zi) > 8A, return to Step g] with b:=2b. 

7. Partition £i , . . . ,£„ into k clusters Ci , . . . ,Q. Compute (well separated) enclosing disks Di , . . . ,Dii 
for the clusters. If the clustering fails to find k clusters and corresponding disks, return to Step|4]with 
b:^2b. 

8. For each /, let A, denote the disk with the same center as D, but with an «-times larger radius. We 
now verify the existence of |C, | roots (counted with multiplicity) of p in A,. If the verification fails, 
return to Step |4] with b := 2b. 

9. If the verification succeeds, output the disks A, (in Step|7] we guarantee that the disks A, are pairwise 
disjoint) and report the number |C, | of root approximations z G {zi , . . . ,z„} contained in the disks as 
the corresponding multiplicities. 

Note that Steps |5]and|6]ensure that logM(z,) = 0{rp + \ogn) for all /, and logH^Li M{zi) = 0{\og{\\p\\ /\pn\))^ 
0(log« + Tp) = 0{nrp). The following Lemma guarantees that the algorithm passes these steps if ^ > b^. 



(9) 

(10) 
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Lemma 5. For any b> bo, it holds that \zi\ <2^^^ for all i, andYYi^yM{zi) < 8A. 

Proof. In the proof of Lemma|4] we have already shown that | < (1 + 1/ (2«^)) •M(z,) for all /. Hence, 
it follows that < (1 + 1 /(2«2)) . 2^i' <2-2^p< 2^p+\ and 

A , X A , 4Mea(n) 4|lp|U 4|lp|l 

il T=l \Pn\ \Pn\ \Pn\ 

□ 

2.2.2 Clustering 

After candidate approximations zi , . . . , z„ are computed using a fixed precision parameter b, we perform a 
partitioning of these approximations into k clusters Ci , . . . , Q, where k is given as an input. The clustering 
is described in detail below. It works in phases. At the beginning of a phase, it chooses an unclustered 
approximation and uses it as the seed for the cluster formed in this phase. Ideally, each of the clusters 
corresponds to a distinct root of p. The clustering algorithm satisfies the following properties; (1) For 
b < b(j, the algorithm may or may not succeed in finding k clusters. (2) For b > bo, the clustering always 
succeeds. 

Whenever the clustering succeeds, the cluster C, with seed f, is contained in the disk D, := A(z,, r,), 
where r,- w min(^, 2^7)' ^^'^ — Furthermore, for b > bo, D, contains the rootz,- (under 

suitable numbering) and exactly m, many approximations. 

For the clustering, we exploit the fact that an approximate factorization p{x) = PnYlii-"^ ~ ^i) of p with 
Hp — pII < 2^^ 11/7 1 1 yields m,- approximations z of the root Zi with distance less than 2^^^^^"''^ to Zi (at least, 
for b > bo). Thus, in simple terms, we aim to determine clusters C of maximal size such that the pairwise 
distance between two elements in the same cluster is less than 2 • 2^''/(^l''l'. We next give the details. 

1. Initialize to the empty set (of clusters). 

2. Initialize C to the set of all unclustered approximations and choose z G C arbitrarily. Let a:— 2 L'°S"J +2 
and 5 := 2-''/\ 

3. Update C to the set of points q gC satisfying \z — q\<2 "^-ifB. 

4. If |C| > a/2, add C to ^€ . Otherwise, set a ; = fl/2 and continue with step[3j 

5. If there are still unclustered approximations, continue with step|2l 

6. If the number of clusters in is different from k, report failure, double b and go back to the factor- 
ization step. 

Note that, for b > bo, the disks A(z, ,2^^/'^'"'') are disjoint. Let Z, denote the set of root approximations 
in A(z,, 2^*''(^'"''). Then, |Z,| = m, according to Lemma 2] We show that, for b > bo, the clustering 
algorithm terminates with C = Zi if called with an approximation zG Zi. 

Lemma 6. Assume b > bo, Zi G Z,-, Zj G Zj, and i ^ j. Then, 

\zi-Zj \ > 2 ■ (2-V(i6'«,) 4_2"''/('^'"J-^). 

Proof Since b > bo, we have 2-^/(2»'f) < for i = i, j by ^ and 2-^/(1^'"') = (2-V(2'nf))i/8 < c7£/(4«) < 
(Jf/8 by (H. Thus, 

\zi - Zj I > max((J„ Oj) - 2-^1 ^^"'-^ - 2-''/(2"';) 

> (Ji/2 + Oj/2 - Oi/A - Oj/A 

> 2 ■ (2"^/('^'"') -(-2"''/(^^"'^'). 

□ 
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Lemma 7. If b > bo, the clustering algorithm computes the correct clustering, that is, it produces clus- 
ters Ci to such that Ci ~ Zi for all i (under suitable numbering). Let Zi be the seed of Q and let 
Gi = m\\\j^i \zi — Zj\. Then, (1 — l/«)c7,- < a, < (1 + 1 /«)(?,■ and Ci as well as the root Zi is contained in 

Proof. Assume that the algorithm has already produced Z\ to Z,_i and is now run with a seed z G Z,-. We 
prove that it terminates with C = Z;. Let ^ be a power of two such that I < m,- < 2£. The proof that the 
algorithm terminates with C = Z,- consists of two parts. We first assume that steps 2 and 3 are executed for 
a = 2£. We show that the algorithm will then terminate with C = Zj. In the second part of the proof, we 
show that the algorithm does not terminate as long as a > 2£. 

Assume the algorithm reaches steps 2 and 3 with a/2 = £, i.e. a/2< nii < a. For any approximation q g 
Zj, we have \z — q\ < 2 • 2^''/^^"''-' = 2 "''^^/5 < 2 "^^[8. Thus, Z\ C C. Conversely, consider any approximation 
q ^ Zi. Then, |z — <?| > 2 • 2^*/''^'"'' > 2 '"'^ifB > 2 and thus no such approximation is contained in C. 
This shows that C = Zi. Since |C| > a/2, the algorithm terminates and returns Z,-. 

It is left to argue that the algorithm does not tenninate before a/2 = £. Since £ and a are powers of 
two, assume we terminate with a/2> 2£, and let C be the cluster returned. Then, m, < a/2<\C\ < a and 
Z, is a proper subset of C. Consider any approximation q £C\ Zi, say q £ Zj with j ^ i. Since q ^ Z,, we 
have \q-z\ > 2 • (2-*/(i^'"') + 2-''/('*^'"-')) > 2 • 2-^/(1^'"') > 2 And since ^ e C, we have \q-z\< 

2 "^yfS. Thus, 4m, < a/2 and, hence, there are at least 3fl/8 many approximations in C\Zi. Furthermore, 

\zi-Zj\ < \zi-z\ + \z-q\ + \q-Zj\ < 2-''/(2'«'' + 2 °^^ + 2-^/(2'"y) < 2-^(16'",) ^2 '"^ + 2-''/('*'«v) < 

3 "'v^. Consequently, there are at least 3fl/8 roots Zj ^ Zi counted with multiplicity within distance 3 "'v^ 
to Zi. This observation allows us to upper bound the value of namely 

iPil^YlUi-Zj^J < (3 ''/^)3"/82(«-"'/-3«/8)r;, 

Mi 

a contradiction to (fTol i. 

We now come to the claims about a, and the disks defined in terms of it. The relation between c, and 
a,- follows from the second part of Lemma|4] All points in C, = Z,- have distance at most 2 • 2^''/'^^'"'^ from 
Zi. Also, by © and Q, 

2 • 2-''/(2m,0 < jnin(l /„2^ (7;/(512«2)) < min(l /n^, ai/{256n^)) 
Hence, C,- as well as Zi is contained in A(z;,min(l/n2, a,/ (256n^))). □ 
Lemma 8. For a fixed precision b, computing a complete clustering needs 0{nb + n^Tp) bit operations. 

Proof. For each approximation, we examine the number of distance computations we need to perform. 
Recall that b (property (|5)) and a are powers of two, a < An by definition, and b>%n> 2a by porperty (|5)- 
Then, "'^f5 = 2-^/(2«) g 2"^. Thus, the number "'^ has a very simple format in binary notation. There 
is a single one, and this one is b/ {2a) positions after the binary point. In addition, all approximations z 
have absolute value less than 2 • 2F due to Step|5]in the overall algorithm. Thus, each evaluation of the form 
|z — ^1 < 2 "^-sfB can be done with 

0{Y + log = 0{{b/a) + r) = 0{{b/a) + Y,, + logn) 

bit operations. 

For a fixed seed z, in the /-th iteration of step 2, we have at most a < n/2'^^ many unclustered approx- 
imations left in C, since otherwise we would have terminated in an earlier iteration. Hence, we perform 
at most a evaluations of the form |£ — ^| < 2 "^\f5, resulting in an overall number of bit operations of 
a ■ 0{{b/a) + r) = 0{b + aF) for a fixed iteration. As we halve a in each iteration, we have at most 
log n + 2 iterations for a fixed z, leading to a bit complexity of (9 (/? log « + nP) = 0{b + nr) = 0{b + nYp). 

In total, performing a complete clustering has a bit complexity of at most 0{nb + n^Tp). □ 
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When the clustering succeeds, we have A: clusters Ci toQ and corresponding seeds Zi , . . . ,Zjt C . . . ,£„}. 
For / = 1, ... we define D, := A(|,,r,), where f, is the seed for the cluster C, and 

r, :=min(2-r2i°g"1^2r'°g'^'/(256«')l) >min(^,^^). (11) 

In particular, r,- is a 2-approximation of min(l/n^, a',/(256n^)). Notice that the cost for computing the 
separations a,- is bounded by d{nb + n^Fp) bit operations since we can compute the nearest neighbor graph 
of the points z, (and thus the values a,) in (9(«log«) steps ifTOl with a precision of 0{b + riF). 

Now, suppose that b > bo. Then, according to Lemma |7] the cluster C, is contained in the disk D,. 
Furthermore, D, contains exactly one root z, of p (under suitable numbering of the roots), and it holds 
that nii = mult(z,,/?) = |C,| and min( 1/(2«2),C7,/(512«2)) < n < min(l/«2,c7,/(64n2)). If the clustering 
succeeds for a b < bo, we have no guarantees (actually, the termination condition in step 4 gives some 
guarantee, however, we have chosen not to exploit it). Hence, before we proceed, we verify that each disk 
Di actually contains the cluster C, . If this is not the case, then we report a failure, return to the factorization 
step with b = 2b, and compute a new corresponding clustering. 

In the next and final step, we aim to show that each of the enlarged disks A,- := A{zi,Ri) A{zi,nri), 
i = l,...,k, contains exactly one root Zi of p, and that the number of elements in C, C A, equals the 
multiplicity of Zi- Notice that, from the definition of r, and A,, it obvious that the disks A, are pairwise 
disjoint and that C, C D, C A,. 

2.2.3 Certification 

In order to show that A, contains exactly one root of p with multiplicity |C, |, we show that each A, contains 
the same number of roots of p and p counted with multiplicity. For the latter, we compute a lower bound 
for \p{z)\ on the boundary bd A, of A,, and check whether this bound is larger than \{p~ p){z)\ for all points 
z G bdA,. If this is the case, then we are done according to Rouche's theorem. Otherwise, we start over the 
factorization algorithm with b = 2b. We now come to the details: 

1. Let X = 2'^ be the 2-approximation of \\p\\ /\p„\ as defined in step|2]of the overall algorithm. 

2. For / = 1,. . . let z* -r,- e A,-. Note that |z*| < (1 + 1/n) ■M(z;) since nn < l/n. 

3. We try to establish the inequality 

\p{z*)/p„\>Ei ■.= 32- 2-' XM{zi)" (12) 

for all i. We will see in the proof of Lemma [TO] that this implies that each disk A, contains exactly 
one root Zi of p and that its multiplicity equals the number |C, | of approximations within A, . In order 
to establish the inequality, we consider p = 1,2,4,8, . . . and compute \p{z*)/pn\ to an absolute error 
less than 2^P. If, for all p < b, we fail to show that \p{z*)/pn\ > Ei, we report a failure and go 
back to the factorization algorithm with b = 2b. Otherwise, let p,- be the smallest p for which we are 
successful. 

4. If, at any stage of the algorithm, p,- > b, we also report a failure and go back to the factorization 
algorithm with b = 2b. Lemma |9] then shows that, for fixed b, the number of bit operations that are 
used for all evaluations is bounded by d{nb + n^Zp + rr'). 

5. If we can verify that \p{zi +nri)/pn\ > Ei for all /, we return the disks A, and the multiplicities 
nii = \Ci\. 

Lemma 9. For any i, we can compute \p{z*i)/ Pn\ to an absolute error less than 2^P with a number of bit 
operations less than 

d{n{n + p+n\ogM{zi) + Tp)). 
For a fixed b, the total cost for all evaluations in the above certification step is bounded by 0{nb + n^Tp + 
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Proof. Consider an arbitrary subset S C {zi, . . . ,£„}. We first derive an upper bound for HrGS k; ~ ^1- 
For that, consider the polynomial ps{x) := HiGsl-^^^z)- The ;-th coefficient of ps is bounded by (''J') • 
X\zesM(z) < 2" Yl"^iM{zi) < 8A • 2" due to step|6]in the overall algorithm. It follows that 

Yl\z*~z\^\ps{z*)\<in + l)M{z*)"-U-2" 

zeS 

<64(n+l)2.2"2^"M(z;)" 

In order to evaluate \p{z*)/pn\ =11^=1 |z* we use approximate interval evaluation with an absolute 
precision K = 1,2,4,8, . . .. More precisely, we compute the distance of z* to each of the points zj, j = 
up to an absolute error of 2^^, and then take the product over all distances using a fixed point 
precision of K bits after the binary point@ We stop when the resulting interval has size less than 2^ P. 
The above consideration shows that all intermediate results have at most 0{n + Tp + nlogM(2,) bits before 
the binary point. Thus, we eventually succeed for an K = 0{p + Tp + n + «logM(z,)). Since we have 
to perform n subtractions and n multiplications, the cost is bounded by 0{nK) bit operations for each K. 
Hence, the bound for the evaluation of \p{z*)/pn \ follows. 

We now come to the second claim. Since we double p in each iteration and consider at most \ogb 
iterations, the cost for the evaluation of \p{z*)/pn\ are bounded by d{n{n + Pi + n\ogM{zi) + Xp)). Since 
we ensure that p,- < b, it follows that the total cost is bounded by d{nb + n'^Xp + rr' + log (nf= i M(z,' ) ) ) . 
The last summandis smaller than «^ • 8A according to step|6] and A < 2 < 2(« + 1)2^''. This shows 

the claim. □ 

We now prove correctness of the certification algorithm. In particular, we show that the inequality (fT2t 
implies that the disk A, contains the same number of roots of the polynomials p and p. 

Leiiuna 10. 1. For all points x e bdA;, it holds that 

w)\>\m)\- 

2. If inequality ( 1721 ) holds for all i, then A,- isolates a root of Zi of p of multiplicity nii = mult(zi,p) = |C,|. 
3- Ifb> bo, then 

|p(zni /min(256,a,.)y-;^^^^^^_,^^^.^„ 



\p„\ V 1024n 

Proof. First, we show that, for any two points x,y e bd A,, their distance to any approximation z G Dj differ 
only by a factor of (1 + 1/n): 

{l-l/n)\y-z\<\x~z\<{l + l/n)\y-z\. 
For a fixed z, assume x to be the farthest point on bdA from £, and let y be the nearest. For ; ^ j, we have 

\x-z\ < \x^Zi\ + \zi-Zj\ + \zj~z\ < + ,and 
\y-z\ > \z-zj\ - \y~zi\ - \zj-z\ > (1 - • 



Similarly, for / = j: 



\x-z\ < \x-Zi\ + \zi-z\ < (1 + l/n)nr; ,and 
\y-z\ > \y-Zi\ - \zi-z\ > (1 - l/n)nn ■ 



Consequently, for any x,y G bdA,, it holds 



\pi^)\ = \Pn\-U\^~U>{l + l/nT\p„\-fl\y-ze\>\\p{y)\ 

e=i (=1 ^ 



*In fact, we compute an interval Ij of size less than 2 ^ such that \z* —Zj \ & Ij, and then consider the product Ylj^j- 
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This shows the first claim. 

We turn to the second claim. Since «r; < 1 /«, we have |x| < (1 + 1 /n)M{zi) for all x e bdA;. Now, if 
\p{.z*)/pn\>2>2-2-''XM{zi)",\htn 

W)\ > \piz*m > 8|p„|A2-^(l - l/«)"M(x)" 

> Ml-'Mix)" > \\P'P\\M{xr > \p{x)~p{x)\. 

Hence, according to Rouche's theorem. A, contains the same number (namely, |C, |) of roots of p and p. If 
this holds for all disks A,-, then each of the disks must contain exactly one root since p has k distinct roots. 
In addition, the multiplicity of each root equals the number |C, | of approximations within A,. 

It remains to show the third claim. Since b > bo, it follows that min(l/(2«^), C7,/(512«^)) < r; < 
min(l/«^, (J,/ (64«^)) and |2, — z, | < r,; cf. the remark following the definition of r, in ( fTTT i. Thus, 

\piz*)\>\p{zr)\^2-''\\p\\-Mizr)" 

= b(z,- + (z,-z,+nr,-))|-2-^||p||-M(zn" 

>{{n-l)nr\p„Pi\/4-4-2-''\\p\\M{zir 
'(n— l)min(256,(7; 
512n 



> ( ^ ^^:rr ) b„P,|/4-4.2-''||p||M(z,r 



where the first inequality is due to \{p ^ p){x)\ < 2 •M(x)", the second inequality follows from 

\zi—Zi + nri\ < (« + l)r,' < a,/«, Lemma[3]and M(7*) < (1 + 1/n) •M(z,), and the third inequality follows 



from r,- > min(i, ^^Si^). In addition, we have 



smce 



1/ II \~,) _ ^Q24„ J 4096 



MMiziT 

< 2-''l^ ■ 2-^1^ ■ 2^1' • (« + 1) •M(z;)" 

/min(256,c7,)\'"' . , x , ^„ 

^ / min(256,g,) Y"' ^ ^ ^min(256, CJ,) ^^ '"' |p„P,| 



V 1024« ) 27«-i - V 1024« ) 4096 



where the second inequality follows from (fTOl i. (|9), and (|6]l Q and the third inequality follows from Tp < 
nVp + n + 1 (LemmalB and M(z,)" < 2""^/'. Finally, 

\piA)\ ^ /min(256,a,.)y'^^^^^^^^_,MMfer 



|;5„| V 1024« ; 8 - \pn\ 
>64-2-^AM(z,)", 

where the first and the second inequality follow from (ITJt and the third inequality holds since A is a 2- 
approximation of \\p\\ / |p„| and |z,f' < (1 + 1 /«)" < 4 □ 

Lemma 11. There exists a b* >bo bounded by 

O (nlog n + nYp + ^J^j (logM(P,- ' ) + logM((J,-' )) 



'Observe 2"*/^"'') < minf-^ , < min(256.cTi) 
VJ05.BIVC z ^ '^2^' Toli?'' - 1024/1 



12 



such that the certification step succeeds for any b > b* . The total cost in the certification algorithm (i.e. for 
all iterations until we eventually succeed) is bounded by 

6 (n^ + n^Tp + n • 1 (logM(Pr ' ) + \ogM{ar ' )) ) 

bit operations. 

Proof. Due to LemmafTOl 

Thus, in order to verify inequality (fTZt . it suffices to evaluate \p{z*i)/pn \ to an error of less than \p{z*) /2p,j\. 
It follows that we succeed for some p, with 

Pi = 0{mi log n + nii max ( 1 , log a;r ' ) + log max( 1 , 1 " ' ) ) . 

In Step 3 of the certification algorithm, we require that the sum over all p, does not exceed b. Hence, we 
eventually succeed in verifying the inequality ( fT2] i for all / if b is larger than some b* with 

b* =0{bo + £ . Mi log n + £ .(logM(P,- 1 ) + m, logM(a;r 1 ) ) ) 
= 0{n \ogn + nFp + £ . {\ogM{Pr ' ) + m,- \ogM{ar ' ) ) ) . 

For the bound for the overall cost, we remark that, for each b, the certification algorithm needs (5(«^ + nb + 
n^Tp) bit operations due to Lemma|9] Thus, the above bound follows from the fact that that we double b in 
each step and that the certification algorithm succeeds under guarantee for all b > b*. □ 

2.3 Complexity of Root Isolation 

We now turn to the complexity analysis of the root isolation algorithm. In the first step, we provide a 
bound for general polynomials p with real coefficients. In the second step, we give a simplified bound for 
the special case, where p has integer coefficients. We also give bounds for the number of bit operations 
that is needed to refine the isolating discs to a size less than 2^'^, with K an arbitrary positive integer 

Theorem 3. Let pix) £ C[x] be a polynomial as defined in Section |Z7] We assume that the number k of 
distinct roots of p is given. Then, for all i = I,. . . the algorithm from Section |Z21 returns an isolating 
disk A{zi,Ri) for the root n together with the corresponding multiplicity mu and Ri < 
For that, it uses a number of bit operations bounded by 

O [n^+nhp + n-Y,l^ {logM{Pr') +milogM{ar'))^ (14) 

The algorithm needs an approximation of precision L of p, with L bounded by 

O (nrp + £j^j (logM(P,-i) +,„,logM((J,-^))) . (15) 

Proof. For a fixed b, let us consider the cost for each of the steps in the algorithm: 

• Steps 1-3, 5 and 6 do not use more than d{n^Tp + nb) bit operations, 

• Step 4 and 7 do not use more than Oln^Tp + nb) bit operations (Corollary [Hand Lemma|8}, and 

• Step 8 and 9 use a number of bit operations bounded by (fl4t (LemmafTTT). 

In addition, for a fixed b, the oracle must provide an approximation of precision (9 (nFp + b) of p in 
order to compute the bound F for Fp, to compute the 2-approximation A of \\p\\ /\pn\, and to run Pan's 
algorithm. The algorithm succeeds in computing isolating disks if b> b* with a b* as in Lemma[TT] Since 
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we double b in each step, we need at most log/?* iterations and the total cost for each iteration is bounded 
by ( fT4l ). This shows the complexity result. 

It remains to prove the bound for When the clustering succeeds, it returns disks Dj ~ A(z,, r,) 
with minfT^T, ^^^-7) < n < mm(-^,-^-^) for all ; = 1, . . . ,m. It follows that Ri = n ■ r, < -rS-, and thus 

\z,-Zj\ > \z,-Zj\ - \zi-z,\ - \zj-Zj\ > \z,~Zj\ • (1 - l/(64«)) > \zi-Zj\/2 for all ij with / ^ j. We 
conclude that a,- > d',72 > 64«/?,. □ 

We remark that the bound (fl4] i can also be reformulated in terms of values that exclusively depend on 
the degree « and the geometry of the roots (i.e. their absolute values and their distances to each other). 
Namely, according to Lemma [T] we have Tp < « + 1 + log ^^j^j^, and the latter expression only involves 
the degree and the absolute values of the roots of p. This yields the bound (|2]i from the introduction. 

In the next step, we show that combining our algorithm with Pan's factorization algorithm also yields a 
very efficient method to further refine the isolating disks. 

Theorem 4. Let pix) be a polynomial as in Theorem\3\ and K be a given positive integer. We can compute 
isolating disks A,(z,-,/?;) with radius Ri < 2^'^ in a number of bit operations bounded by 

3§ + d{nK- maxmi), (16) 

l<z<A: 

where S§ is bounded by ( 1741 ). For that, we need an approximation of precision Lof p with L = Ji' + 0{nK- 
maxi<,<jt'«i)' where J£ is bounded by lilSi . 



Proof. As a first step, we use the algorithm from Section l272] to compute isolating disks A,- = A(z,,/?,) with 
Ri < <7,/ (64n). Each disk A, contains the root z,, nii = mult{zi,p) approximations z e {£1, ...,£„} of z,, and 
it holds that (7,/2 < a, < 2 a,. Let 

We claim that 1 /2|P,| < < 2\Pi\. Since \zi - Zi\ < (7,/(64«) for all /, it holds that (1 - ^^)\zi - z/| < 
l^; — ■?! < + '^,)\zi ~ Zj \ for all j 7^ / and £ G A,-. Thus, |P/ 1 is a 2-approximation of \Pi\. Similar as in the 
certification step, we now use approximate interval arithmetic to compute a 2-approximation jij of and 
thus a 4-approximation of A completely similar argument as in the proofs of Lemma|9]and LemmafTTI 
then shows that we can compute such /x,'s with less than 0{n^ +n^Tp + «LilogM(Pr')) bit operations. 
Now, from the 2- and 4-approximations of C7, and we can determine a b^ such that 

• the properties (|5]l to (|8) are fulfilled, and 



Then, from Corollary [T]and Lemma|4] we conclude that Pan's factorization algorithm (if run with b > b^) 
returns, for all /, m; approximations z of Zi with |£ — z,| < 2^''/(^'"'' < 2^'^. Thus, for each ;, we can simply 
choose an arbitrary approximation z G A, and return the disk A(z,2^'^) which isolates z,-. The total cost 
splits into the cost for the initial root isolation and the cost for running Pan's Algorithm with b = bK- Since 
the latter cost is bounded by d{nbf: + n^Tp), the bound ( fTSI l follows. □ 

Finally, we apply the above results to the important special case, where p is a polynomial with integer 
coefficients. 

Theorem 5. Let p{x) e Z[x] be a polynomial of degree n with integer coefficients of size less than 2^. Then, 
we can compute isolating disks A{zi,Ri), with Ri < for all roots zt together with the corresponding 
multiplicities mi using 

d{n^+n^T) (17) 

bit operations. For a given positive integer K, we can further refine the disks A; to a size of less than 2^^ 
with a number of bit operations bounded by 

d{n^+n^T + nK). (18) 
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Proof. In a first step, we compute the square-free part p* = p/ gcd{p,p') of p. According to 1201 §1 1.2], we 
need d{n^x) bit operations for this step, and p* has integer coefficients of bitsize 0{n + t). The degree of p* 
yields the number A; of distinct roots of p. Thus, we can directly apply our algorithm from Section lZ2l to the 
polynomial p. In order to derive the bound in (ITTt . we have to reformulate the bound from ( fT4l i in terms of 
the degree n and the bitsize T of p. From Jg] Theorem 2] , we conclude that m, log max ( 1 , C7,^ ' ) = <5 («^ + 
«t). Furthermore, we have Tp < T. Hence, it remains to show that ^^^^ logM(Pj^') = (5(n^ +«^t). For that, 
we consider a square-free factorization p{x) ~ YYi=\ [Qti^))' with square-free polynomials Qi G such 
that Qi and p/Q\ are coprime for alH = 1, . . . ,«. Note that the roots of Qi are exactly the roots of p with 
multiplicity I, and that Qi is a constant for most /. We further denote p :— p/\cf{p) and Q[ :~ G//lcf(2/). 
Let Si denote the set of roots of Qi. Then, from the definition of P,, 

n i^'i- n Y{\^i-'^r 

= u u n n 1^'--^/ 

^Y\\{p/Q\){zi)\ - n 

= |res(/5/e|,e,)|.|res(a,(a)')|' 
|res(jV6;,6/)| 

|lcf(e/)"-'*ga(lcf(;,/e|))degaf 

res(e/,e',) 



> 



Icf (2/ )2 deg a - 1 (deg e, )deg 2/ 
1 

lcf(e/)|""'-|lcf(/7)|'l'=g2; -n'dega 



where ms{f,g) denotes the resultanj^ of two polynomials / and g. For the last inequality, we used that 
i'ss{p/Q'i,Qi) € Z and res (2/, 2/) <= Z- Taking the product over all / yields 



ni^-i> 

i=l 

> 



U1=i(icf (Qi)"-' ■\cf{p)<^^sQi -n'-^^sQi) 



|lcf(p)p"-«" 



>2 



-2«T — ttloSH 



Note that, for any /, we also have 



niilpn 



mi\\p„\ 



and thus, 

i k 
^logM(P,.-') =0(«T + n^logM(z,)) =(5(«t), 

where we used that L(logM(z,)) < logMea(/5) < log \\p\\ < log(« + 1) + T. This shows ( fTTI i. 

For the bound in ( fTSl l for the cost of refining the isolating disks A, (z, ,/?,) to a size of less than 2^'^, 
we consider the square-free part p* . Note that the disks A, obtained in the first step are obviously also 
isolating for the roots of p* {p and p* have exactly the same distinct roots) and that Rj < a{zi,p) /(64«) = 



For univariate polynomials 



ref (/, = M{ff-^^ M{g)'^^f H - >') 

(.v,y):/(jr)=«(>')=0 

.v:/(;i)=0 
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a{zi,p*)/ {64n) < a{zi, p*)/ {64degp*). Thus, proceeding in completely analogous manner as in the proof 
of Theorem|4](with the square-free part p* instead of p) shows that we need 0{rr' +rp-X + nK) bit operations 
for the refinement. This proves the second claim. □ 



3 Curve Analysis 

In this section, we show how to integrate our approach to certify roots of a univariate polynomial in an 
algorithm to compute a cylindrical algebraic decomposition JTKTTKlJllISlClEllHlIlll- More specifically, we 
apply the results from the previous section to a recent algorithm, denoted TOPNT, from f3l for computing 
the topology of a real planar algebraic curve. This yields bounds on the expected number of bit operations 
for 

• computing the topology of a real planar algebraic curve, and 

• isolating the real solutions of a bivariate polynomial system 

which improve the currently best bounds ||8] [14] from 0{n^T + n^T^) (deterministic) to 0{n^ + n^z) (ran- 
domized) for topology computation and from 0{n^ +rPz) (deterministic) to 0{n^ +n^T) (randomized) for 
solving bivariate systems. 

3.1 Review of the Algorithm TopNT 

The input of the algorithm is a bivariate polynomial / e 'Z[x,y] of total degree n with integer coefficients 
of magnitude bounded by 2^. The polynomial defines an algebraic curve 

C:^{{x,y)eC^:f{x,y)^0}CC\ 

The algorithm returns a planar straight-line graph embedded in that is /ioto/^/tQ to the real part 
CM:=CnR2 ofC. 

In the first step (the shearing step), we choose an i e Z at randorrPl and consider the sheared curve 

C, := {[x,y) e C2 : ^x.y) := f{x + s-y,y) = 0}. 

Then, any planar graph isotopic to the real part C, r Cs n R of Cs is also isotopic to Cr, and vice versa. 
We choose s such that the leading coefficient (with respect to y) of the defining polynomial fs{x,y) of Cs 
is a constant. This guarantees that Cj r has no vertical asymptote and that it contains no vertical line. By 
abuse of notation, we write C ~ Cs and / = /j throughout the following considerations. 

In the projection step, the x-critical points of C (i.e. all points (xcjo) G C with fyixQ.yo) = 0, fy := ^) 
are projected onto the real x-axis by means of a resultant computation. More precisely, we compute 

. 7?:=res(/,/,;y)eZ[x], 

• its square-free part R* := R/ gcd{R,R'), 

• isolating intervals /i , . . . ,7^ for the real roots ai , . . . , of /?*, 

• the multiplicity m,- := mult(a,-,7?) of a, as a root of R for all ; = 1, . . . ,m, and 

• arbitrary separating values j3o, . . . , A„+i £ K with < j3,„+i, and j3,_i < a, < j8, for all / = 1, . . . ,m. 
We further compute 

• A* — gcd(i./v) fy gcd{i,/,)' 

'We actually consider the stronger notion of an ambient isotopy, but omit the "ambient". is ambient isotopic to Cm if there is a 
continuous mapping : [0, 1] x n- with 0(0, ■) = id]j2, 0(1,Ch) = and 0(<o,-) is a homeomorphism for each /q £ [0, 1]. 
'"initially, we consider s = 0. 
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• G:=res(/;,/^;3'), and 

• the multiplicity Z, := mult(a,,2) of a, as a root of Q for all / = 1, . . . ,m. 

In the lifting step, we compute the fibers of C at the points a, and j3„ that is, we isolate the roots of the 
polynomials faXy) '■— f{oLiiy) S Rb] and f^.(y) f{Pi,y) G ^[y]- For that, we first compute the number 
of distinct complex roots of each of the latter polynomials, and then use the root isolator from Section |2] 

Obviously, each polynomial fp. (y) has A:(j3,) = degfp. = n distinct complex roots. The difficult part is 
to determine the number k{a) of distinct roots of fa{y) for a root a of R* . According to 13] (3.6)] and 13] 
Theorem 5], we have 

k+{a) :=«-mult(a,/?)+muh(a,e) <k{a), (19) 

and, for a generic shearing factor s (i.e. for all but n'^^'' many s), the equality A:+(a) = k{a) holds for all 
roots a of R* . Summation over all complex roots of R* then yields 

K+ -.^ k+{a)^n-dtgR* -dtgR + Aeggcd{R°°,Q) 

a:R'(a)=0 

> Y ^i^) —'■ ^tid K = for generic 

where gcd(/?°°, 2) is defined as the product of all common factors of R and Q with multiplicities according 
to their occurrence in Q. The crucial idea is now to compare the upper bound with a lower bound K 
which also equals AT up to a non-generic choice of some parameters. In order to understand the computation 
ofK^, we first consider the exact computation of K: Let Sres, (/,/y;3') <E "^[x^y] denote the i-th subresultant 
polynomial of f and (with respect to y), and sr,(x) := sres,(/,/y;3') e Z[x] its leading coefficient. In 
particular, we have R = sreso{f, g;y) = res(/,/,,;y). We define; 

So~R*, 5,:=gcd(5,_i,sr,) (20) 

^ _ ^0 _ So ^ _ 5,_i _ gcd(5o,...,5,_i) 
Si gcd(5o,5i)' Si gcd{So,...,Si-i,Si)' 

where ; = !,...,«. Then, ni>i^i constitutes a factorization of R* such that /?,(a) =0 if and only if 
f(ci,y) has exactly « — / distinct complex roots; see 13] Section 3.2.2] for details. Hence, we have K = 
Li>i ~ ■ degRi. Unfortunately, this computation of K is costly in practice due to the computation of 
the full subresultant sequence. Instead, in order to derive a lower bound , we do not carry out the above 
computations directly in Z but in a modular prime field. More precisely, we choose a prime p at random, 
compute the modular images sr|''^(ji:) = STeSi{f mod p,fy mod p;y) S 'Zp[x\ of STi{x) GZ[x], and perform all 
computations from ( |20] i in Zp[x]. This yields polynomials R^-''^ e Zp[x]. Now, |j3] Lemma 4] shows that 

/:-:=^(n-/).deg/?<''^</i:, (21) 

and K = K for all but finitely many bad primes^ Hence, if K < , we have either chosen a bad prime 
or a bad shearing factor In this case, we start over with a new s and choose a new prime p in the lifting 
step. If K = K^, we know for sure that = K, and thus A:+(a) = k{a) for all roots a of the resultant 
polynomial 7?. 

We can now use our method from Section |2]to isolate all complex roots of the fiber polynomials /«, (y) 
and fp. {y). Namely, we can ask for arbitrary good approximations of a, and j3,- (by refining corresponding 
isolating intervals), and thus for arbitrary good approximations of the coefficients of the fiber polynomials. 
In addition, we know the exact number of distinct roots of either polynomial. From the isolating regions 
in C, we then derive isolating intervals for the real roots together with corresponding multiplicities. If one 
of the polynomials /«. (y) has more than one multiple real root, we start over and choose a new shearing 
factor s. Otherwise, we proceed with the final step. 



"in the computation of the 5,'s and R,'s (over Z), all intermediate results have integer coefficients of bitsize bounded by (jit)"'''. 
Since the product of W distinct primes is larger than N\ = 2'^('^'°s'v)^ there exist at most (ht)"''' many bad primes for which ^ K. 
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Connection step. We remark that, except for finitely many s, each has exactly one multiple root. 
The previous two steps already yield the vertices of the graph Namely, these are exactly the pointQ 

y(^) := {(x,y) G : 3/ with x = Ui orx = /3;, and f{x,y) = 0} 

Since each polynomial fa {y) has exactly one multiple root, there exists a unique vertex v along each vertical 
line, where either the number of edges connecting v to the left or to the right may differ from one. Hence, 
connecting all vertices in an appropriate manner is straightforward; see ||3] Section 3.2.3] for more details. 



3.2 Complexity Analysis 

Throughout the following considerations, we say that a multivariate polynomial G E , . . . with 
integer coefficients has magnitude (A^, ) if the total degree of G is bounded by and all coefficients have 
absolute value of 2^* or less. In addition, we fix the following notations: For an arbitrary a G C, 

• we define fa{y) YJi=ofa.iy' '■— f{(^,y) G C[y], where / G '^[x,y] is our input polynomial. We fur- 
ther define log max, \fa^i | > as the logarithm of the maximal absolute value of the coefficients 
of/a- 

• the number of distinct roots of /„ is denoted by k{a). We further denote Za.i 7 ■ ■ ■ i^a k(a) distinct 
roots of fa, and niaj, i = 1, . . .,k{a), the corresponding multiplicities. 

• Ga.i is the separation of Za.i, and Paj = UjM^"'' ~ ■ 

• For an arbitrary polynomial G G C[.y], we denote V{G) the set of all distinct complex roots of G, 
and ^(G) the multiset of all complex roots (i.e. each root occurs a number of times according to its 
multiplicity). 

We first prove the following basic result: 

Lemma 12. For a fixed positive integer k, /ef G G , . . . be an integer polynomial of magnitude 
(N,^). Then, each divisor g G , . . . o/G has coefficients ofbitsize 0{pi +N). 

Proof. We prove the claim via induction over k. For a univariate G G we remark that Mea{g) < 

Mea(G) < ||Gj|2 < 2^^^N, and thus the absolute value of each coefficient of g is bounded by 2^Mea(g) < 

For the general case, we write 

g{xi,...,Xk) = ax{xk)x\^ ■■■xl'i^,W\thaxE'L[xk]. 

For a fixed G {0, . . . ,A^}, the polynomial § (jc 1, . . . G Z[xi , . . . ,JCjt_i] is a divisor of G{xi,. . . ,xi^_i,Xk) G 

Z[jt:i, . . . Since \x^\'^ <N^ = 2^^^°^^ andax{xk) has degree or less, it follows that G(xi , . . . 

has bitsize 0{NlogN + ji). Hence, from the induction hypothesis, we conclude that the polynomial 
g{x\,. . . ,Xk-i,xi;) has coefficients of bitsize (5(/i + A^), and thus ai{i) G Z has bitsize (5(/i + A^) for all 
i = 0,... ,N and all A. Since ai is a polynomial of degree at most A^, it follows that a^ is uniquely deter- 
mined by the values ax{i), and thus Lagrange interpolation yields 

^ x-(x-l)---{x-i+l){x-i-l)---{x^N) 

Expanding the numerator of the fraction yields a polynomial with coefficients of absolute value 2'^^^^°s'^\ 
and thus each coefficient of a^ (xj-) has bitsize d{lJ.+N) because ai (/) has bitsize d{lJ.+N) and there are 
N +1 summands. This proves the claim. □ 



'^For a graph with rational vertices, you may replace each xq = (Xj (or xq = Pi) by an arbitrary rational value in its corresponding 
isolating interval, and the same for each real root of (y). 
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We now come to the complexity analysis for TopNT. For the shearing step, we remark that there exist 
at most n"^''' many bad shearing factors s for which our algorithm does not succeed; see JS] Thm. 5] 
and 111 Prop. 1 1.23]. Thus, when choosing s at random, we can assume that we succeed for an integer s of 
bitsize 0(log«). It follows that the sheared polynomial /(jc + sv,}') has magnitude («, (9(T + log«)). Hence, 
throughout the following considerations, we can assume that the leading coefficient of / (with respect to y) 
is an integer constant and that / has magnitude («, (9(t + log«)). 

Lemma 13. We can compute the entire subresultant sequence Sres,(/,/,,;y), with / = 0, . . . ,n, the poly- 
nomial Q ~ Tes{f* ,f*;y), and the square-free parts R* and Q* of the corresponding polynomials R = 
Sieso{f,fy;y) —ies{f,fy;y) and Q with (5(«^ + «^t) bit operations. 

Proof. For two bivariate polynomials g,h e Z[x,y] of magnitude (A^,/x), computing the subresultant se- 
quence SKSi{g,h;y) £ Z[x,y] together with the corresponding cofactor representations (i.e. the polynomi- 
als Ui,\'i € with Ujg + Vjh = Sresi{g,h;y)) needs d{N^jJ.) bit operations ifTTl l5l. The total degree 
of the polynomials Sres, (/,/,, is bounded by A^^, the y-degree is bounded by A^— /, and all coefficients 
have bitsize 0{Nli). Furthermore, according to ll20] §1 1.2], computing the square-free part of a univariate 
polynomial of magnitude {N,jJ.) uses d{N^jJ.) bit operations, and the coefficients of the square-free part 
have bitsize 0{N + li). Hence, the claim concerning the computation of the polynomials Sres, (/,/j,;y) and 
R* follows from the fact that / and fy have magnitude («, (9(t + log«)) and R has magnitude (n^, (5(«t)). 

From ||2] Prop. 10. 14, Cor.10.15], we conclude that the polynomials f* and /* can directly be obtained 
as cofactors in the subresultant sequence of f^ and /,,. Thus, their computation needs 0{n^z) bit operations 
as well. Since /* divides and /* divides /,,, Lemma[T2lyields that f* and /* have magnitude {n,0{n + 
t)). Thus, for computing Q, we need d{n^x) bit operations, and Q has magnitude {n^,d{n^ +«t)). Its 
square-free part Q* can be computed using bit operations 0{n^ i + n^). □ 

We now bound the cost for computing and comparing the roots of R and Q. 

Lemma 14. The roots of the polynomials R and Q can be computed with d{n^ + n^z) bit operations. The 
same bound also applies to the number of bit operations that are needed to compute the multiplicities 
mult{a,R) and mult(a, Q), where a is a root ofR. 

Proof. According to Theorem |5] we can compute isolating disks for the roots of the polynomials R and Q 
together with the corresponding multiplicities with d{n^ + n^x) bit operations since R and Q have magni- 
tude {n^,0{n^ +nT)). For each root a of R, the algorithm returns a disk A'*'(a) := A{a,ra) with radius 
ra < 64 "e^^ ; and thus we can distinguish between real and non-real roots. A corresponding result also 
holds for each root j3 of Q, that is, each j3 is isolated by a disk A^^^fi) with radius less than Fur- 
thermore, for any given positive integer K, we can further refine all isolating disks to a size of less than 2^*^ 
with 0{n^ + n^T + n^K) bit operations. 

For computing the multiplicities mult(a, Q), where a is a root of R, we have to determine the common 
roots of and Q. This can be achieved as follows: We first computed := deggcd{R* ,Q*) for which we need 
d{n^ +«^t) bit operations. Namely, computing the gcd of two integer polynomials of magnitude (A^,jLi) 
needs 0{N'-pL) bit operations. We conclude that R and Q have exactly d distinct roots in common. Hence, 
in the next step, we refine the isolating disks for 7? and Q until there are exactly d pairs (A(*'(a),A(2)(j3)) 
of isolating disks that overlap. Since P . = R ■ Q has magnitude {2n^,0{n^ +nz)), the minimal distance 
between two distinct roots a and j3 is bounded by the separation of P, thus it is bounded by 2^'-"'" +" 
We conclude that it suffices to refine the isolating disks to a size of 2^*^'" +"^), hence the cost for the 
refinement is again bounded by (5(«* + «^ ) . Now, for each of the d pairs (A('^)(a),A(2)(j3)) of overlapping 
disks, we must have a ^ j5, and these are exactly the common roots of R and Q. □ 

From the above Lemma, we conclude that we can compute the numbers (a) for all roots a of R with 
0{n^ +n^z) bit operations. Thus, the same bounds also applies to the computation of the upper bound 

~ Y,a for K = ^(oc)Ill For the computation of the lower bound , we have the following 

result: 

"For simplicity, we ignored tliat (in practice) can be computed much faster from tlie equality = n degR* — degR + 
deggcd{S°°,g) instead of computing tlie first and, tlien, summing up all values. 
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Lemma 15. Wfe can compute K with d{ir'x) bit operations. 

Proof. We have already computed the leading coefficients sr, e of the subresultant sequence Sres, (/,/,, ;y) 
and the square-free part Sq := R* of the resultant polynomial R. Note that all polynomials 5, and Rj as 
defined in ( |20l ) have coefficients of bitsize d{nT + n") because all of them divide R*. Thus, except for 
0{nT) many bad primes, the modular computation over Zp yields polynomials ^I'^^T?!'^' G with 

deg7?|'^' = degRj for all /, and thus = K. Hence, we can assume that we only have to consider primes p 
of bitsize (9(log(nT)). Since we can compute the polynomials sr,- and R* using (5(«^t) bit operations, the 
same bound also applies to their modular computation over ZnF^ 

For the computation of the polynomials 5 ■ G Zp [x] , we have to perform at most n gcd computations 
(over Zp) involving polynomials of degree n^. Thus, the cost for these computations is bounded by 0{n ■ 
n^ log (« t) ) bit operations since computing the gcd of two polynomials in Zp [x] of degree can be achieved 

with d{N) arithmetic operations in Zp due to ll20l Prop. 11.6]. For the computation of the 7?|'^''s, we 
have to consider the cost for at most n polynomial divisions. Again, for the latter computations, we need 
0{n ■ n^log(«T)) bit operations. □ 

We remark that it is even possible to compute K directly in an expected number of bit operations 
bounded by O(n^T). Namely, computing the gcd of two integer polynomials of magnitude {N,li) needs 
an expected number of bit operations bounded by d{N^-\-Nli) according to ||20l Prop. 11.11]. Hence, this 
yields the bound (5(«(«^ + n^x)) for the expected number of bit operations to compute the polynomials Si. 
Obviously, the same bound also applies to the computation of the 7?,'s. 

For the analysis of the curve topology algorithm, it remains to bound the cost for isolating the roots 
of the "fiber polynomials" fa/iy) G M.[x] and fp.{y) G M.[x], where the a,'s are the real roots of R and the 
j3,'s are arbitrary separating values in between. In practice, we recommend to choose arbitrary rational 
values )3,, however, following this straight forward approach yields a bit complexity of (5(n^ + n^z) for 
isolating the roots of the polynomials fp. (y) G Q[y]. Namely, if j3, is a rational value with of L,-, then fp. has 
bitsize d{nLj + t). Thus, isolating the roots of fp. needs 0[n^Li + n'^x) bit operations. However, since the 

separations of the a,'s are lower bounded by 2^'^("*+"'^\ we cannot get anything better than d{n'^ + n^x) 
for the largest L,-. 

The crucial idea to improve upon the latter approach is to consider, for the values j3,, real roots of the 
polynomial R{x) instead, where R := gcd{{R*y (r*)") defined as the square-free part of the derivative of 
R*. Note that the polynomials R and R do not share a common root, and, from the mean value theorem, we 
further conclude that, for any two consecutive real roots of R, there exists a root of R in between these two 
roots. We can obtain such separating roots by computing isolating disks for all complex roots of R such 
that none of these disks intersects any of the isolating disks for the roots of R. The computation of R needs 
0{n^ +n^T) bit operations since (R*)' has magnitude {n^,0{n^ +nz)). We can use the same argument as 
in the proof of Lemma[T4]to show that it suffices to compute isolating disks for R and R of size 2^'^^" +" 
in order to guarantee that the disks do not overlap. Again, Theorem |4] shows that we achieve this with 
0{n^ + n^z) bit operations. 

Now, throughout the following considerations, we assume that the separating elements j3, are real roots 
of R with j8,_ 1 < a, < Pi. We will show in Lemma[T9lthat. for isolating the roots of all polynomials fp. and 
/a,, we need only 0{n^ + n^z) bit operations. For this purpose, we need the following result: 

Lemma 16. Let G G Z[x] be a polynomial of magnitude (N,^). For an arbitrary subset V' C Y{G), it 
holds that 

^ logMea(/„) = diNi + npL), and ^ Ta = (9(A^T + njU). 

aeV aeV 

In particular, for G G {R,R}, the bound write as (9(n^ +«^t). 



'''Note that, in practice, we never compute tlie entire subresultant sequence over Z. In the proof of Lemma[T3] we only assumed 
their exact computation in order to keep the argument simple and because of the fact that our overall complexity bound is not affected. 
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Proof. The proof is almost identical to the proof of Lemma 5 in lfT4l . The only difference is that we 
consider a general G, whereas in |fT4), only the case G = R has been treated. Note that Meaa > 1 for 
every a G V{G), and that the Mahler measure is multiplicative, that means, Mea(g)Mea(/i) = Mea{gh) for 
arbitrary univariate polynomials g and h. Therefore, 

5^ logMea(/„) < logMea(/„)=logMea( !]/«)■ 

aeV aei'(G) \aer(G] ) 

Considering / as a polynomial in x with coefficients in 1\y\ yields 

n /« = ^T^^7^^ ^ logMea(/„)<logMea(res(/,G;x)). 

It is left to bound the degree and the bitsize of res(/, G;x). Considering the Sylvester matrix of / and G 
(whose determinant defines res(/, G;x)), we observe that it has n rows with coefficients of G (which are 
integers of size 0{iJ,)) and rows with coefficients of / (which are univariate polynomials of magnitude 
(«, T + log«)). Therefore, the y-degree of Tes{f,G;y) is bounded by 0{nN), and its bitsize is bounded by 
+log«) +A^(T + log«)) = d{Nz + niJ.). This shows that logMea(res(/,G;x)) = d{Nz + nli), and 
thus the first claim follows. 

For the second claim, note that the absolute value of each coefficient of fa (y) is bounded by (« + !)• 
XM{a)'\ where A = 2'^(^+>°g") is an upper bound for the absolute values of the coefficients of /. Thus, we 
have 

Lt„< Y L log((« + l)AM(a)«) 

aeV aer(G) aeflG) 

= C»(A^(T + logn) +«logMea(G) = d[N'Z + n^) 
For the last claim, note that, for G G {R,R\, we have < «^ and n = (5(«^ + nt). □ 
Lemma 17. For G e {R,R\, we have 

k(a) 

Y L ma,,logM(cT„ j) = d{n' + nh) . 
aeV{G) i=i 

£ £logM(P„j) = 0(n4 + „\). 

aeV{G) i=i 

Proof. First, consider G =R. For any root a of R, we define m{a) := mult(a,/?). From (fT9] l, we conclude 
that L/i'i'('na,/ —!)=«- k{a) < « — k{a) +mult(g, a) < mnlt{a,R). Furthermore, since m{a) > 1, it 
foUows that niaj < m{a) + 1 < 2m{a) for all /. Hence, we get 

k{a) k{a) 

Y i:m„,,logM(cT„j)< Y 2.m(a)-^logM(cJ„j) 

aeV{R) 1=1 aeV(R) i=l 

^^ol Y ^a)-(«logMea(/„)+logM(sr„_,(„)(a)-i)) ) 
io( nlogMea(/„)+logM(sr„_,(„)(a)-i) ) 

For (1), we used lfl4l 9] to show that 

m 

Y logM((J„j) =0(«logMea(/„)+logM(sr„_,(„)(a)-i)) . (22) 

1=1 
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(2) follows from the fact that each a G V{R) occurs m{a) times in y{R). Finally, for (3), we apply 
Lemma[T6lto bound the first sum and lfT4l Lemma 8] to bound the second one. 
The second claim can be shown as follows. For each a, we first spht the sum 

k{a) k{a) 

£logM(P„j)= |^log|P„,r'+ I iog\Paj\. (23) 

'■=1 '=1 i;\Paj\>l 

Then, for the first sum, we have 

k{a) k{a) 

^ log|P„,r' - L \ogY\\Za.-Za.j\-""^'' 

i=l i=\ j^i 

/k(a) k(a) 

< log ™^ 1 rXpi + y riogcT^h^-'' 

(2) _ 

= (9(nT)+log|sr„_^(„)(a)| ' + ^ ^oga^\Y,{ma.i - I) 

'=1 

(3) '^(«) 

< (5(«T)+log|sr„_^„)|-i +m(a)- ^ logM((7„j) 

(=1 

^= (5(«t) + (m(a) + 1) • (nlogMea/a +logM(sr„_i,(„))-') (24) 

For (1), we have rewritten the product as a subresultant term, where we used |l2] Prop. 4.28]. Furthermore, 
the distances \za.i — Za.j \ have been lower bounded by the separation of Za,i- For (2), note that Icfa is an 
integer of bitsize (9(T + log«), that k < n, and that Wiiiiai < «". For (3), we used that Ly!=i'' ~ 1) !i 
m{a), and, in (4), we applied (|22] |. Now, summing up the expression in (l24l) over all a G V(^) yields 

^ ';£\og|p„,r'-(5(«4+«'T), 

aGV(«) (=1 

where we again use Lemma[T6land ||l4l Lemma 8]. 

For the second sum in ( |23] |. we use that (cf. proof of Theorem|5) 

ma,i\\cf[fa) 



and thus 



/ / k(a) 

52 E logl^«,'i='5 E «Ta+«logn^(Za,<) 

= (5( E («T„+«logMea(/„)) ) 



according to Lemma [17] We conclude that 



E 'Eiog^(^ai)-^(«'+«'^)- 

aeV{R) 1=1 
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Now, we consider the case G = R. Note that, for each a E V{R), we have niaj = 1 for all /, and thus 
k{a) = n. Namely, R and R do not share a common root, and thus each polynomial fa has only simple 
roots. Also, V{R) ~ 'f{R) since R is square-free. The following computation now shows the first claim 

k(a) 

L £ma,,logM(cT„j) 

aeV(R) 1=1 

= (hogM{a-,])\ 

aev{R) \'=l / 

=(9 «logMea(/a)+logM(sro(a)-') 

\aeV(R) 

aeV{R) 

In order to bound the sum in the above expression, note that 

£ logM(/;(a)-')= Y \og\R{a)\-'+ ^ \R{a)\. 
aeV{R) aeV{R) a;\R{a)\>l 

We first compute an upper bound for each value |/?(oc)|. Since R{x) has magnitude (n^,(5(«T)), it follows 
that \R{(x)\ has absolute value less than 2'^("^' ■M{a)" . Hence, for any subset V' C V{R), it follows that 

Y \og\R{a)\ < d{n\) +n^logMeaiR) = Oin'^ + n\). 

aeV 

Thus, it is left to show that La log l-'^lof)! ^^ = 0{n'^ + n^t), which follows from 

£ log|i?(a)r'=log n \R{a)\-' 
aeV(R) a£V(R) 

\_ \res{R,R)\ J 

In the second equation we rewrote the product in terms of the resultant res{R,R) ||2l Prop. 4.16]. Since 
R and R have no common root, we have |res(/?,^)| > 1. Thus, the last equation follows from the fact that 
the leading coefficient of R has bitsize (5(«^ + nz) and that deg{R) < n^. 

Similarly, for the second claim, we first derive an upper bound for Y.aev{R) I^(:fa;>i log l^a,ii- Again, 
we can use exactly the same argument as for the case G = R to show that the latter sum is bounded by 
d{n'^ + n^z). Hence, it suffices to prove that 

aeV(R) '=1 
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This result follows from 



k(a) 

^ i:iog|p„,r' 



aeV(R)i=^ J^i 




|sro(a)| 

|lcf(/a)2"-2| 
|lcf(/a 



|2n-2| 



The last step follows from dZST l. We remark that above computation is similar to the one for the case G=R. 
However, we used the fact that R is square-free, and thus all multiplicities ma,i are equal to one. □ 

Lemma 18. Let G G {R,R} and let La (zN be arbitrary positive integers, where a runs over all real roots 
ofG. Then, we can compute approximations of precision La for all polynomials f{cc,y) using 



din^ + nh + n^Y.^a) 



bit operations. 



Proof. For each a, we use approximate interval arithmetic to compute an approximation of the polyno- 
mial fa. If we choose a fixed point precision p, and a starting interval of size that contains a, then 
the so-obtained interval approximation of fa has interval coefficients of size 2-P+2(« + l)22W(a)"; see 
again fiy., Section 4] and [|T4l Section 5] for more details. Thus, in order to get an approximation of preci- 
sion La of fa, it suffices to consider a p of size (5(T + «logM(a) +La). Thus, by doubling the precision 
p in each step, we eventually succeed for some p = pa = (5(T + «logM(a) +La). The cost for the inter- 
val evaluations is then dominated (up to a logarithmic factor) by the cost in the last iteration. Thus, for a 
certain a, the cost is bounded by (5(«^(t + «logM(a) +La)) since, for each of the « + 1 coefficients of 
fa, we have to (approximately) evaluate an integer polynomial (i.e. the coefficients of / considered as a 
polynomial in y) of magnitude {n,0{T + log«)) at x = a. The total cost for all a is then bounded by 

oin^- T + «logM(a)+L„ I =(5(n*'+n5T + «'^La), 

where we again used the result in LemmafTT] For the interval evaluations, we need an approximation of the 
root a to an absolute error of less than l^P". Such approximations are provided if we compute isolating 
disks of size less than 2^^ for all roots of G, given that jcis larger than max^ Pa = d{T + nmaxa logM(a) + 
maxaLa) = d{n^ + n^T + maxaLa). In the proof of Lemma [T4l we have already shown that we can 
compute such disks using 0{n^ + n^T + n^K) bit operations. Thus, the claim follows. □ 

Lemma 19. Let G € {R,R}. Then, computing isolating disks for all roots of all fa, (X G PR, together 
with the corresponding multiplicities uses 

d{n^+nh) 

bit operations. 

Proof. For a fixed a, let Ba be the number of bit operations that are needed to compute isolating disks for 
all roots of fa together with the corresponding multiplicities. According to Theorem|3] we have 

Ba = («3 + nha + n ■ J^Jf"' (m„,, logM((J„ j ) + logM(P„ | ; 
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The corresponding algorithm from Section |Z2] returns isolating disks for the roots Za,i and their muti- 
plicities niaj- Furthermore, since the radius of the disk isolating za,i is smaller than Ga.i/ (64«), we can 
distinguish between real and non-real roots. The algorithm needs an approximation of precision La of fa 
with 

La = d («T„ + {ma, logM((T„ ) + logM(P„ J )) ) . 

From LemmafTSl we conclude that we can compute corresponding approximations for all /«, OC G V{G) n 
R, with a number of bit operations bounded by 

0(«'' + «^T + «2. £ La). 

aeV(G)(M 

The above expression is bounded by d{n'^ + n^x) because 

E («^«+LJ^i' ('«a,,logM(cT„j)+logM(P„j))) 

aeV(G) 

is bounded by 0(n^ + tt't) according to Lemma [16] and [17] The same argument also shows that the sum 
over all Ba is even bounded by Oiii^ + ri^x). Hence, the claim follows. □ 

We can now formulate our main theorem: 

Theorem 6. Computing the topology of a real planar algebraic curve C — V{f), where f e 1j\x^y\ is a 
bivariate polynomial of total degree of n with integer coefficients of magnitude bounded by Y , needs an 
expected number of bit operations bounded by 

d{n^ + nh). 

Proof. We already derived a bound of 0{n^ + n^x) or better for each of the steps in the projection and in 
the lifting phase of our algorithm. The final connection phase is purely combinatorial since we ensure that 
each fa, with a a root of the resultant R, has at most one multiple real root. Thus, we can compute all 
adjacencies in linear time with respect to the number of roots of critical and intermediate fiber polynomials. 
Since their number is bounded by 0{n^), this step can be done in 0{n^) operations. □ 

Finally, obtain the following result for solving bivariate polynomial systems: 

Theorem 7. Let g, h G ^[x,y] be coprime polynomials of magnitude («, x). Then, we can compute isolating 
boxes for the real solutions of the system g{x,y) — h{x,y) = with an expected number of bit operations 
bounded by 

d{n^ + nh). 

Proof. The idea is to consider the polynomial f{x,y) := g^ + h^ and to compute the topology of the curve 
C := C]R defined by /. Since g and h are assumed to be coprime, the system g = h = has only finitely 
many solutions, and the set of these points coincides with the "curve" C. Hence, the topology algorithm 
returns a graph that consists of vertices only. According to Theorem|6] the cost the topology computation 
is bounded by (5(n^ + n^x) bit operations in expectation since / has magnitude {2n,0{x)). 

However, in general, our algorithm does not directly return the solutions of the initial system but the 
solutions of a sheared system g{x + sy,y) = h{x + sy,y) ~ 0. Here, i is a positive integer of bitsize (9(log«) 
for which TOPNT succeeds in computing the topology of the sheared curve C := Cjjj defined by f{x,y) ~ 
f{x+sy,y) ~ 0. Since C consists of isolated singular points only and there are no two co vertical points (note 
that our algorithm only succeeds for an s for which there are no two co vertical extremal points), it follows 
that, for each point {x,y) € C, i is a root of the resultant R = res(/,/;y) and y is the unique (multiple) 
real root of f{x,y). The point {x,y) is represented by an isolating box B{x,y) = I{x) x I{y), where I{x) 
is the isolating interval for the root x of ^ and I{y) is the isolating interval for the root y of Each 
solution (x,y) of the initial system can now be recovered from a unique solution {x,y) E C. More precisely, 
x=x — s-y and y = y. However, in order to obtain isolating boxes for the solutions {x,y), we have to refine 
the boxes B(x,y) first such that the sheared boxes B{x,y) := (/(i) —s ■ I{y),I{y)) do not overlap. Note 
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that the latter is guaranteed if both intervals I{x) and I{y) have width less than cj(x,^)/(4|s|) < a{x,R)/4. 
Namely, if the latter inequality holds, then the intervals I{x)—s- I{y) are pairwise disjoint. Hence, it follows 
that the corresponding isolating intervals have to be refined to a width less than w{x,y) = (7{x,R)/n'^^^\ 
For the resultant polynomial R, we conclude from Theorem |5] that computing isolating intervals of size 
less w{x,y) uses 0{n^ +n^z) bit operations since \ogM{w{x,y)^^) — (5(«'*+«^t) and ^ has magnitude 
{n^ ,0{nz)). In order to compute an isolating interval of size w{x,y) or less for the root y of f{x,y) (in fact, 
for all roots of f{x,y)), we need 

d{n' + nhs + n- ^'^"J (m^,, logM(c7r/ ) + logM(P^.i )) 
+ («maxm;i- ,) • logM(w(i,3')^') 

bit operations; cf. the proof of Lemma [19] with a = x and f ^ f- Also, we need an approximation of 
precision Lf of f{x,y) with Lf bounded by 

0{nT, + Y!iS {mAogMiaiJ ) + \ogM{Pi/ )) 
+ (nmaxm^,,) • logM(w(i,3')^'). 

Since max, m^,,- < 2 •mult(i,^) and w{x,y) = a{x,R)/n^^^\ it holds («max,mi,,) •logM(w(i,y)^') = d{n ■ 
mult(i,^) • logM(c7(i,^)^'). Thus, summing up the cost for computing the roots of f{x,y) over all real 
roots of R yields the bound d{rr' +n^x). Here, we use an analogous argument as in the proof of Lemma[T9] 
and the fact that £f « • mult(i,^) • \ogM{a{x,R)^^ = 0{n^ + n^x). The more costly part is to compute the 
approximations of precision Lf of the polynomials Again, we can use Lemma [T6l and [TTl to show 

that Lr = (5(n* + tv't). Thus, from Lemma [Ts] we conclude that the approximations of the /(i,y)'s can 
be computed with d{n^ + rr'x) bit operations. □ 
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